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ABSTRACT 



O: 
U 

I I We use XMM-Newton observations of the galaxy cluster 

O , MaxBCG J217. 95869+13. 53470 to analyze its physical properties and dynami- 

cal state. MaxBCG J217.95869+13.53470 is found at a redshift of 0.16, has a mass of 
^1 ~ 1 X 10^* Mq, and a luminosity of 7.9 x 10^'^ ergs~^ . The temperature map shows 

the presence of hot regions towards the north and west of the brightest cluster galaxy 
(BCG). From the entropy distribution, regions of high entropy match the location 
! of the hot regions; more high entropy regions are found to the west, and ~ 165 kpc 

■ to the southwest of the central AGN. A second X-ray bright galaxy is visible ~ 90 

' kpc to the northeast of the BCG, at a redshift of 0.162. This galaxy is likely to be 

the BCG of a smaller, infalling galaxy cluster. The mass of the smaller cluster is 
~ 10% the mass of MaxBCG J217. 95869-1- 13.53470, yielding an impact parameter 
of ~ 30 — 100 kpc. We compare the results of our X-ray observations with GMRT 
. observations of the radio source VLSS J1431. 8-1-1331, located at the center of the 

' cluster. Two sources are visible in the radio: a central elongated source that bends 

at its northern and southern ends, and a southwestern source that coincides with a 
region of high entropy. The radio sources are connected by a bridge of faint radio 
emission. We speculate that the southwestern radio source is a radio relic produced 
by compression of old radio plasma by a merger shock. 



00 

o 



X 



\-{ . Key words: galaxies: clusters: individual: MaxBCG J217. 95869+13. 53470 - X-rays: 

galaxies: clusters ~ shock waves 



1 INTRODUCTION 

Studies of large-scale structure (LSS) formation have found 
that nearly all massive clusters have undergone several merg- 
ers in their history, and that present clusters are still in the 
process of accreting matter. A significant fraction of the ac- 
creting matter is in the form of smaller clusters and galaxy 
groups. 

Observations have shown that 10 — 35% of massive 
clusters host diffus e radio sources (|Giovannini et al] 1 19991 : 
IVenturi et all |2007| . [ 20081 ). In general, these diffuse radio 
sources are only found in clusters that display significant 
substructure in the gala xy distribution and in t he thermal X- 
ray emitting gas (e.g., ISchuecker et al.ll200ll : IVenturi et al.l 
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|2008| ). These clusters are also characterized by high X-ray 
luminosities and temperatures. All of these signatures are 
indicative of the cluster undergoing a merger. 

The diffuse radio sources observed are also referred to 
as radio halos or relics. Radio halos result from the reac- 
celeration of a relatively old electron population as a re- 
sult of galaxy cluster mergers. Radio relics are irregularly 
shaped radio sources with sizes ranging fr om 50 kpc to 2 
Mpc. They can be divided into two groups (|Kempner et al.l 
|2004| ). Radio gischt are large elongated, often Mpc-sized 
radio sources located in the periphery of merging clus- 
ters. They probably trace shock fronts in which parti- 
cles are accelera te d via diffusive shock accele r ation (PSA ; 
e.g., lBel]|[T978a l'bl: 'Bla ndford fc O strikci^'l9 Z^; lDrurvlll98l : 
iBlandford fc Ei chlcr 1983; IJones fc Ellison Il99j ). Accord- 
ing to DSA, if the particle acceleration is still ongoing, the 
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integrated radio spectrum should follow a single power-law. 
The second category of relics are radio phoenices. In these 
objects, fossil radio plasma from a former activity cycle of a 
FRI/II radio galaxy is compress e d by a merger shock wav e 
lEnfilin fc Gopal-Krishnal I2OOII : lEnfflin fc Briiggenl l2002h . 
When a fossil radio cocoon is passed by a cluster merger 
shock wave, the cocoon is compressed adiabatically and not 
shocked because of the much higher sound speed within it. 
Therefore, shock acceleration cannot be the mechanism that 
re-energizes the relativistic electron population. However, 
the energy gained during the adiabatic compression com- 
bined with the increase in the magnetic field strength can 
boost the radio emission from the fossil radio cocoon. One 
prerequisite for this is that the electron population is not 
older than 0.2 — 2 Gyr, depending on the amb ient pressure 
of the cocoon jEnfilin fc Gopal-Krishnal I2OOII '). The radio 
plasma in such a radio phoenix has a steeply curved spec- 
trum due to synchrotron and inverse Compton ( IC) losses. 
Propo sed e xamples of these are th ose found by ISlee et al] 
(I2OOII ') and (|van Weeren et al.ll20lil ). 

The radio source VLSS J1431.8-f 1331, located in the 
galaxy cluster MaxBCG J217.95869+13. 53470, at a redshift 
of 0.16, is a good candidate for a shocked radio plasma 
due to its extreme spectral index and disrupted morphol- 
ogy. VLSS J1431. 8-1-1331 was included in a sample of dif- 
fuse ste ep spectrum sources sele cted from the 74 MHz VLSS 
survey (|van Weeren et alll2009l ). The spectral index of the 
source is very steep with a = —2.03 ± 0.05 between 74 and 
1400 MHz. GMRT observations showed the source to consist 
of two main components (see Fig. [l}. The brightest eastern 
component is associated with the BCG and can be classi- 
fied as a radio AGN. For the western component no obvi- 
ous optical counterpart is detected. The western component 
is located at a distance of about 165 kpc from the clus- 
ter centre and has a largest linear size (LLS) of 125 kpc. 
The radio plasma from the AGN appears to be affected by 
the ICM because it bends eastwards at the northern and 
southern ends of the sou rce. GMRT 325 MHz observations 
(Ivan Weeren et aLll201ll ) reveal a faint bridge of radio emis- 
sion connecting the two components suggesting a link be- 
tween them. The spectral index of the radio core of the 
AGN is —0.5. The spectral index steepens to a < —2 at the 
northern and southern ends of the radio AGN. The western 
component has a very curved radio spectrum with a — —1.5 
between 325 and 610 MHz, steepening to a = —2.5 between 
610 and 1425 MHz. The very steep curved radio spectrum 
indicates that the radio plasma has undergone significant 
synchrotron and inverse Compton losses. It is likely that 
this radio plasma originated from the central AGN given 
the faint radio bridge connecting the two components. 

Here we report results of XMM-Newton observations 
of MaxBCG J217.95869-H3.53470, obtained with the Euro- 
pean Photon Imaging Cameras (EPIC). The spectra allow 
for a detailed spatial analysis of the distribution of physi- 
cal properties (e.g. temperature, electron number density, 
pressure, entropy) within the central region of MaxBCG 
J217.95869-f 13.53470. Comparing the X-ray results with 
the results obtained from GMRT observations of the ra- 
dio source allows us to study the interaction between the 
BCG, the radio plasma and the ICM, as well as the cluster's 
merging history. 

Section[2]describes the observations and the X-ray data 



reduction process, while Section[3]presents the results of the 
spectral analysis. Section |4] provides an interpretation of the 
X-ray and radio results. Finally, Section [S] summarizes the 
paper. 

Unless stated otherwise, we use a flat ACDM universe 
with = 0.3, Q.A = 0.7, and Ho = 70 kms"^ Mpc"^ . At 
a redshift of 0.16, 1 arcsec corresponds to 2.758 kpc. 



2 DATA REDUCTION AND SPECTRAL 
ANALYSIS 

2.1 X-ray data reduction 

On 2009 May 25, MaxBCG J217.95869+13.53470 was ob- 
served for 45 kiloseconds (ks) with the EPIC MOS and pn 
cameras of XMM-Newton. The detectors were operated in 
full-frame mode, using the medium filter. 

The data was reduced using the XMM-Newton Science 
Analysis System (SAS), version 9.0.cQ. For each of the three 
detectors, we extracted two lightcurves: one in the hard 
band, 10 - 12 keV (12 - 14 keV for pn), and one in the 
soft band, 0.3 — 10 keV. These light curves were binned to 
100 and 10 seconds, respectively. Time periods when the 
count rate deviated from the mean by more than Sa were 
excluded from the data. The higher bin size used for binning 
the hard band light curves is due to the lower count rate in 
this energy band. Nevertheless, excluding periods with high 
count rates in the hard band eliminates most of the flaring. 
Repeating the procedure in the 0.3 — 10 keV energy band 
further elimantes any short flares that might have been pre- 
viously missed. The resulting net effective exposure times 
were 41 ks for MOSl, 41.3 ks for MOS2, and 35.7 ks for 
pn. 

Background subtraction was performed using blank sky 
observations developed at Leicester Universi ty by the EPIC 
Blan k Sky team (based on the work of ICarter fc ReadI 
|2003), and fiher- wheel closed (FWC)0 datasets, which were 
filtered in the same way as the target datasets. A composite 
background for each detector was created from the blank 
sky and FWC datasets by normalizing the FWC event rates 
to the difference between the target and blank sky event 
rates outside the field of view in the energy range 10 — 12 
keV (MOS) and 12 - 14 keV (pn). This compensates for dif- 
ferences in the particle background between the target and 
background observations. 

As customary, out-of-time events were scaled by a factor 
of 6.3??lf| and subtracted from the pn data. 

A vignetting-corrected and background-subtracted flux 
map of the cluster, created by combining the images in the 
0.4 — 7 keV band from all three detectors, is shown in Fig. [U 
overlaid with GMRT 325 MHz radio contours. The X-ray 
emission is elongated in the NE-SW direction, parallel to 
the radio emission. The BCG appears slightly offset from the 
radio AGN. The bright source in the northeast is identified 
in the SDSS DR7 as SDSS J143153. 78-^133315.1, an AGN 

^ http:/ /heasarc. gsfc.nasa.gov/docs/xmm/abc/ 

http:/ /xmm2.esac.esa.int/external/xmm_sw_cal/background/ 

filtcr_closed/index.shtml 

http:/ /www. mpc.mpg.de/xray/wavG/xmm/cookbook/ 

EPIC-PN/ootcvonts.php 
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Figure 1. Surface brightness map of MaxBCG 
J217.95869+13. 53470 in the 0.4 - 7 keV energy range. The 
image is eorrected for vignetting and the background was 
subtracted. GMRT 325 MHz radio contours are drawn at 
4xrmsx[l,2,4, 8, ...], with rms = 1.3212 X 10"* Jybeam-^. 



at a redshift of 0.16, the same as the cluster. This AGN was 
excluded from the spectral analysis. 



2.2 Temperature map 

To obtain a more detailed image of the physical proper- 
ties of the cluster, we binned the data using an adaptive 
binning method base d on weighted Voronoi tessellations 
iDiehl fc Statieill2006l). which is a generalization of the algo- 



Cappellari fc CopinI (|200a ). This allows us to create 



rithm of 

maps of the temperature, pressure, and entropy distribu- 
tions, which are presented in the following subsections. 

We binned the combined counts image to a target 
signal-to-noise per bin of 31 (following background subtrac- 
tion) in the energy range 0.4 — 7 keV, which corresponds to 
~ 960 counts per spatial bin. Target and background spec- 
tra, redistribution matrices (RMF) and ancillary response 
files (ARF) were created for each resulting bin. The target 
spectra were grouped to a minimum of 30 counts per spec- 
tral bin. 

Using XSPEC, we fitted the spectra with a Galactically- 
absorbed single-temperature APEC model, which descibes 
emission from coUisionally-ionized diffuse gas. The hydro- 
gen column density, nu, was fixed to 1.52 x 10^" cm~^ - the 
weighted average hydrogen column density at the position 
of the BCG in the Leiden/ Argentine/Bonn (LAB) Survey 
of Galactic Hi. The temperatures and metallicities were left 
free in the fit, but they were coupled to be the same for all 
detectors. The normalizations were allowed to vary indepen- 
dently. The temperature uncertainties in the fitted bins are 
between 8% and 20% (at the 1-sigma level), while the re- 
duced lie between 0.5 and 1.6 (typically higher at larger 
radii, where the emission from larger spatial bins is not well- 
described by a single-temperature model). 

The projected temperature map, with radio contours 



overlaid, is presented in Figure (2] Hot regions are present 
to the north and west of the BCG. It is interesting to note 
that the northern hot plasma is parallel to the northern edge 
of the central radio source, while the western hot plasma is 
approximately coUinear with the south-western radio emis- 
sion. 



2.3 Pressure and entropy maps 

We used the results of the spectral fits to create quasi- 
deprojected pressure and entropy maps of the cluster. 

We calculated the electron number density correspond- 
ing to each bin using the definition of the spectrum normal- 
ization, 



A^bin = 



47r [Da{1 + z)Y 

10-'* 
47r [Da{1 + z)]' 



ricnuVhin , 



(1) 



where Da is the angular size distance, Da = 568.9 Mpc, nn 
is the hydrogen number density, Vbin is the bin volume, and 
we assume rio/riH « 1.2. 

Since the normalization parameters were left to vary 
independently in our previous fits, we refitted the spectra 
with a Galactically-absorbed single-temperature APEC model 
in which the temperature and the metallicity were fixed to 
the previously determined values, while the normalizations 
were left free to fit, but were constrained to be the same for 
all detectors. 

To estimate the volume of each bin, we assume that the 
observed emission originates only from the region between 
the spheres of minimum and maximum radii tangent to the 
bin (|Mahdavi et al.ir2005l 'l. Let -Rmax and -Rmin be the maxi- 
mum and minimum radii, respectively. The volume of a bin 
is then 



V = 



2SL 



(2) 



where S is the projected area of the bin in the sky plane, 
and L — 2^ -Rmax — ^min i^ ^ii^ longest length through the 
contributing volume. The projected area can be easily cal- 
culated by directly counting the number of pixels in the bin, 
and using the fact that the pixel scale is 4 arcsec/pixel. With 
this method of describing the 3D geometry of a bin, we are 
considering a constant temperature along the line of sight. 

The electron number density map is shown in Figure [51 
At the 1-sigma level, the uncertainties are 3-5%. 

Using the bin electron number density and temperature, 
we calculated the pressure, p = n^hTx, and entropy, K = 
kTxric^^^ , where k is Boltzmann's constant, and Tx is the 
temperature of the X-ray emitting gas. The two maps are 
shown in Figure[2] The uncertainties at the 1-sigma level are 
9-20% and 8-20%, respectively, for pressures and entropies. 

A pressure gradient is seen in the pressure map, with 
the pressure decreasing uniformly toward the outskirts. The 
entropy distribution follows a less clear pattern. The south- 
western radio source coincides with a region of higher en- 
tropy. High entropy regions are also present west, northeast 
and northwest of the BCG. 

Figure 13] shows the temperatures and entropies, plotted 
against bin number, in most of the bins covering the center 
of the cluster. The bin numbers are shown in Figure [4] The 
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Figure 2. (a) Temperature map of MaxBCG J217.95869+13. 53470 in units of koV. (b) Electron number density (xlO^^ cm"^). 
(c) Quasi-deprojccted pressure map ( X 10^"'^"'^ Ba). (d) Entropy map ( keV cm^ ). North is up, east to the left. Radio contours as in Figure[l] 
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Figure 3. Left: Temperatures in the bins shown in Figure|4] Right: Entropies in the bins shown in Figure|4] Error bars are plotted at 
the Icr level. 



bins were selected based on the variations observed in the 
temperature and entropy maps. 



3 RESULTS 
3.1 Mass 

To calculate the mass profile of the cluster, we extracted 
spectra from concentric sector annuli centered on the BCG 
(Fig.[5|. This sector was chosen to be separated from the hot 
regions found in the temperature map. As before, the spec- 
tra were grouped to a minimum of 30 counts per spectral 



A merger-revived radio phoenix in MaxBCG J217. 95869+13. 53470 5 




Figure 4. Bin numbers of some of the central regions in the 
cluster. Temperatures and entropies for these bins are shown in 
Figure [3] Radio contours as in Figure [T] 




Figure 5. Concentric sector annuli used for creating temperature 
and electron number density profiles. The sector was chosen such 
that it does not cross the hot regions visible in the temperature 
map of the cluster. The position of the AGN that was not included 
in the analysis is marked by a green circle; the extracted region 
has a radius of 20 arcsec. 



bin, and fitted with an absorbed single-temperature APEC 
model. Tiie results of these fits are summarized in Table [l] 
The no rmalization of the model, together with the volume 
matrix (|Kriss et al.lll98^ allow us to create the deprojected 
electron number density profile of the cluster using the geo- 
metrical deprojection model of lEttori et al.l (j2002h . 



The deprojected electron number density profile was fit- 
ted with a standard beta-model, yielding the best-fit param- 
eters /3 = 0.5 ± 0.02, rc = 88 ± 8 kpc, no = 0.006 ± 0.0004 
cm~^ , while the projected temperature profile was fitted 



Table 1. Radii, temperature, normalization of the best-fit APEC 
model, and reduced of the APEC fit for the sector annuli 
regions shown in Figure [5] Errors quoted at the la level. 



Radius (arcsec) 


T (keV) 


X 10^ 


xVd.o.L 


0-20 


o 9+0.60 
■J-^-0.48 


5.2 ± 0.22 


0.66 


20 - 35 


o ,,+0.79 
■^"^-0.55 


4.8 ±0.23 


0.52 


35 - 55 


9 c+0.39 
^■"-0.30 


8.3 ±0.31 


0.49 


55 - 80 


o 1 +0.51 
■^■'--0.49 


7.5 ± 0.30 


0.48 


80 - 110 


9 c:+0.48 
^■0-0.44 


6.7 ±0.30 


0.64 


110 - 160 


1 o+O.lO 
-'-■'-'-0.10 


12 ± 0.55 


0.80 


160 - 210 


n Qa+0 069 
0.y»_o.070 


8.7 ±0.53 


0.69 


210 - 285 


'J-y^-o.oei 


6.6 ±0.42 


0.87 


with a model of the form: 

rp _ To 



1 + ^ 



(3) 



The best-fit parameters for the temperature model are To = 
3.3 ±0.4 keV, rs = 326 ±196 kpc, 5 = 0.9 ±0.5. The profiles 
and the fitted models are shown in Fig. [6l 

We first determined the mass of the cluster by integrat- 
ing over the electron number densities in the sector in Fig. O 
up to i?5oo = 527 kpc (corresponding to the radius where the 
average density is 500 times the critical density of the uni- 
verse). The gravitational mass is given by M = Mgas/Jb, 
with the mass of the gas Mgas ~ J HemunedV , where 
fie ~ 1.17 is the mean molecular weight per electron, and 
the baryo n fraction U s ing th e gas mass fraction inside 
R500 from I Vikhlinin et al.l (|2006|). and correcting for stellar 
components ( Gonzalez et al.l " 2007^ . the gravitational mass 
of the cluster within -R500 is 9.2 x 10^'' M©. If the baryon 
fraction of the cluster approaches th e mean cosmic baryon 
fraction (0.167; iKomatsu eral]|2010l ') ■ then this sets a lower 
limit of 5.8 x 10^^ Mq on the gravitational mass. 

By assuming hydrostatic equilibrium and spherical sym- 
metry, the gravitational mass of the cluster can also be calcu- 
lated from the models fitted to the electron number density 
and temperature profiles via: 



M(r) 



kTr 
G^mu 



dloe. 



+ 



dlogT 



d log r d log r 



(4) 



where T is the X-ray temperature of the gas at radius r, 
k is Boltzmann's constant, G is the gravitational constant, 
H ~ 0.6 is the mean molecular weight of the ionized plasma, 
and mu is the hydrogen mass. The resulting mass profile is 
shown in Fig. [T] The mass enclosed in a radius of iisoo is 
5.9 x 10" Mq. 

Alternatively, the mass can be estimated from the 



M5, 



Yx scaling relation, where Yx = Afgas,5oo75oo, where 



T500 is the X-ray tempera ture of the gas in the [0.15 — l]f?5oo 
region (|Pratt et al.ll2009l ) of the sector showed in Fig. [S] 

Cluster spectra were extracted for all detectors from 
radius 0.15i?5oo up to i?5oo in the sector in Fig. [5] and re- 
binned to a minimum of 30 counts per spectral bin, to al- 
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Figure 6. Deprojected electron number density (top) and pro- 
jected temperature (bottom) profiles of the cluster, extracted 
from the sector shown in Fig. [5] The models fitted to these profiles 
are shown in red. Error bars arc plotted at the Icr level. 




low for statistics. From the fits to the spectra, which 
were done following the same method as for the bin spec- 
tra, we find a temperature T500 = 2.1 keV, while the de- 
projected electron number density profile yields a gas mass 
A4as,50o = 1.2 X 10^^ Mq within R5QQ. Therefore, using 
equation (1) in iPratt et al.l (|2009l ) . the gravitational mass 
of the cluster within .Rsoo is ~ 1.1 x 10^'' M©. 

Spectra were also extracted from the full sector in Fig. 
[S]up to -R500 (i.e., from radius to -R500), and fitted with an 
absorbed APEC model in the same way as described in Section 
[2:2] The fit yields a fiux of ~ 1.1 x 10"^^ erg cm"^ s"^ in the 
0.5 — 2 keV energy range, equivalent to a luminosity of the 
cluster of 7. 9 x 10"*^ erg s~^ . U sing the L[0.5 - 2] - M500 re- 



rg s . U s: 

lation from lPratt et al.1 (|2009l ). this luminosity corresponds 
to a total mass within i?5oo of 1.8 x 10^"* Mq. 

The disturbed temperature distribution suggests that 
MaxBCG J217.95869-I-13. 53470 is possibly a system un- 
dergoing a merger. By selecting the extraction region of 
our spectra such as not to include the hot northern and 
western regions, our aim was to minimize the influence of 
a possible merger on the different mass estimates. How- 
ever, the assumptions of hydrostatic equilibrium and spher- 
ical symmetry do not necessarily hold. The mass cal- 
culated under the assumption of hydrostatic equilibrium 
could be under estimated due to residual non-thermal sup- 
port (see e.g., IZhang et al]|2008l ). The mass inferred from 
the M500 — Vx and L[0.5 — 2] — M500 scaling relations 
could presumably also be affected by the interaction of 
MaxBCG J217.95869-I-13. 53470 with another galaxy clus- 
ter. Nevertheless, we believe the mass of the cluster to be 
closer to 1 X 10^^ Mq , higher than the mass calculated under 
the assumption of hydrostatic equlibrium. 

3.2 Central cooling time 

The cooling time of optically thin gas is given by 

_ SnkTx , . 

(7-lKnHA(rx,Z) ' 

where n is the total number density (n « 2.3nH), 7 is the 
adiabatic index (7 = 5/3), and A(rx, Z) is the cooling func- 
tion for temperature Tx and metallicity Z. Ado pting the 
cooling function from [Sutherland fc Dopital(|l993h . the cen- 
tral cooling time of the cluster is found to be 6.1 ± 0.53 
Gyr. This indicates that MaxBCG J2 17.95869-H3. 53470 is 
a weak cool core or transition cluster (jHudson et al.|[2010l ). 

However, owing to the point-spread functions (PSFs) of 
the XMM-Newton telescopes, the central X-ray luminosity 
may be more peaked, which would lower the central cooling 
time. Therefore, 6.1 Gyr is only an upper limit. Observations 
with the Chandra telescope would be necessary to see if 
MaxBCG J217.95869-I-13. 53470 is a stronger cool core. 



Figure 7. Gravitational mass profile of the cluster within -R500, 
derived from the fits to the temperature and electron number 
density profiles. 



4 DISCUSSION 

4.1 The dynamical state of the cluster 

To search for evidence of X-ray excess in the central region 
of the cluster, we calculated an azimuthally-averaged surface 
brightness profile by considering 50 8"-wide rings centered 
on the brightest X-ray pixel. This profile was used to con- 
struct a symmetric cluster image, which was subsequently 
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Figure 8. After subtraction of a symmetric brightness profile, 
excess X-ray emission can be easily observed to the NE of the 
cluster center. Overlaid in red are SDSS optical contours. The 
northeastern AGN was excluded from the image, and the excluded 
region is marked by a green circle of radius 20 arcsec. 



subtracted from the original image of the cluster. The result 
of this subtraction is shown in Fig. [H] After subtraction, an 
X-ray excess can be seen ~ 90 kpc to the NE of the cluster 
center. This excess is associated with a few galaxies visible 
in the SDSS optical image, and suggests an interaction of 
MaxBCG J217.95869-I-13. 53470 with another galaxy cluster 
or group of galaxies. 

This hypothesis is further supported by the disturbed 
temperature distribution, and by the bent radio emission, 
all of which point toward a disturbed, merging system. 

To estimate the mass of the second galaxy clus- 
ter, an azimuthally-averaged surface brightness profile of 
MaxBCG J217.95869-f 13.53470 was created while mask- 
ing the emission from the second cluster. Subtracting this 
profile from the original image allows us to put a lower 
limit on the luminosity of the second cluster. The ratio 
of the luminosities of the two clusters is ~ 6. Assuming 
the clusters are in the same dynami cal state an d usin g 
the L500 — M500 scaling relations from iPratt et al. I (|2009D . 
the mass of the second cluster is 10 — 40% of the mass of 
MaxBCG J217.95869-I-13. 53470, with the lower value corre- 
sponding to the scaling relation of a disturbed system. 

Assuming the merging process began at z=0.2 and the 
two clusters dominate the mass in their region of the Uni- 
verse, the greatest separation between them is do ~ 2.5 Mpc. 
The exact age of the Universe at the time of the merger has 
very little effect on the relative velocity an d on the impac t 
parameter. Using equations (23) and (19) in lSarazi 3 (|2002D . 
the impact parameter of the smaller cluster can be estimated 
via 

fc^A^^l^l-A^ "^S{M,,M,), (6) 
where A is the spin parameter (|Peeblesl 19691 ) . d is the current 



separation between clusters, and 
(Ml + Mif 



/(Mi,Af2 



^3/2^^3/2 



(Ml + M2)5/3 



3/2 



(7) 

where Afi and M2 are the masses of the two clusters. 

The relative velo city of the clu sters can be estimated 
from equation (14) of lSarazinI ()2002h : 



2G (Ml + Ma 



1 1 

d do 



1 - 



(8) 



We assume a spin parameter A = 0.05 (|Sarazinll2002l l. 
If the mass of MaxBCG J217.95869-f 13.53470 is ~ 1 x lO" 
M0 and if the current separation between the clusters' cen- 
tres is equal to the projected distance between them, d ~ 90 
kpc, the impact parameter for a disturbed system with a 
mass ratio of 1:10 is ~ 30 kpc, while the relative velocity 
of the clusters is ~ 3200 kms~^ . This velocity is approxi- 
mately six times higher than the sound velocity within R500 
{cs = 568 kms^^ ), and this is possibly only a lower limit 
since the luminosity (hence, the mass) of the smaller clus- 
ter is likely to be underestimated. It is possible, however, to 
overestimate the relative velocity due to projection effects. 
If we assume that the temperature increase observed in the 
north equals the temperature jump across a merger shock, 
then the shock has a Mach number of 1.3, which corresponds 
to a shock velocity of ~ 1000 kms^^ . For this velocity to 
be comparable to the relative velocity of the merging clus- 
ters, the separation between them should be ~ 700 kpc. This 
would yield an impact parameter of ~ 100 kpc. 

We speculate that the smaller cluster entered the ICM 
of MaxBCG J217.95869-hl3.53470 from southeast, giving 
rise to the high-temperature structures visible north and 
west from the BCG. If this is true, and the separation be- 
tween the BCGs of the two clusters is small, then we are 
likely observing the system shortly after closest encounter. 



4.2 Radio relic 

The southwestern radio source is not associated with X-ray 
or optical emission. However, there is a bridge of faint radio 
emission connecting it to the central radio source, suggesting 
that it is related to the radio AGN. The radio emission of the 
southwestern source is elongated and bent eastward, which 
implies that the relic is being acted upon by the surround- 
ing ICM. Its very curved radio spectrum, with a = —1.5 
between 325 and 610 MHz, steepening to a = —2.5 be- 
tween 610 and 1425 MHz, indicates that the plasma has un- 
dergone significant synchroton and inverse Compton loses. 
Unless the source has been "revived" by the passing of a 
shock capable of compressing the faded plasma, the south- 
western source would probably n ot be visible even at the 
relativ ely l ow GMRT frequenc i es. [ Enf^lin fc Gopal-Krishnal 
(I2OOII ') and lEnfilin fc BriiggenI (|2002l ') have shown that the 
energy gained during adiabatic compression, combined with 
the increase in magnetic field strength, can cause the aged 
plasma to emit radio waves again. Given the signs of merg- 
ing present in the X-ray observations, we speculate that 
the southwestern radio source originated as a bubble from 
the central AGN, lost significant energy via synchroton 
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and inverse Compton emission, and was subsequently re- 
vived by a merger shock triggered by the interaction be- 
tween MaxBCG J217.95869+13. 53470 and a smaller cluster 
or group of galaxies. Therefore, we classify the source as a 
radio phoenix. 

The simplest velocity estimate of a rising bubble can be 
obtained by equating the ra m pressure and buoyancy forces 
acting on a bubble (see, e.g., IGuU fc Northoveilll973l ). This 
yields a terminal velocity of the bubble of 

where r is the bubble radius, g the gravitational acceler- 
ation, and C the drag coefficient. The expression for the 
terminal velocity can be further rewritten using the Kep- 
lerian velocity at a given distance from the cluster center: 
V ~ [{r / R){%/'iC)Y^^ vk , where R is the distance from the 
center, and vk = \/gR is the Keplerian velocity. For a solid 
sphere moving through an incompressible fluid, the drag co- 
efficient C is of the ord er of 0.4-0.5 for Reynolds numbers in 
the range of 10^ - 10^ l|Landau fc Lifshitz|[l959l , see, e.g., ). 
Thus, a large and strongly underdense bubble will rise with 
a velocity comparable to the Keplerian velocity. 

In MaxBCG J219. 95869-1-13. 53470 the radio source in 
the SW is at a distance of ~ 1 arcmin from the central 
AGN which corresponds to a distance of 165 kpc. The mass 
enclosed within this radius, calculated by integrating over 
density, is 1.2 x 10^^ Mq, so so the Keplerian velocity 165 
kpc the BCG is vk ~ 550 kms~^ . Therefore, it would take 
at least 300 Myr for the bubble to rise to its current location. 
However, the synchrotron lifetime of the electrons is less 
than approximately 100 Myr - shorter than the buoyant 
rise time. 

In cluster centres, Mach numbers of merger shocks are 
low and hence compression factors inside shock fronts are 
unlikely to exceed 2 — 3. The compression caused by the 
merger shock wave gives rise to a burst of low frequency 
emission, but practically no high frequency emission. This is 
due to the rapid decay of the upper end of the electron spec- 
trum in the dense cluster centres, which essentially wipes out 
the adiabatic energy gains of these electrons. The source de- 
cays on a time-scale of a few tens of Myr, mostly due to the 
heavy synchrotron losses. 

The passage of a merger shock wave through 
the intracluster medium causes irreversible transfor- 
mations that lead to an increase in entropy. In 
MaxBCG J217.95869-f 13. 53470, the entropy map shows an 
overall gradual increase in entropy from the centre toward 
outskirts. However, on top of this natural entropy gradient, 
there are several areas of higher entropy surrounding the 
central radio source. One of these areas coincides with the 
position of the southwestern radio source, and it is possibly 
a result of the processes leading to the revival of old plasma. 



5 CONCLUSIONS 

SDSS optical images of MaxBCG J217.95869-hl3.53470 re- 
veal the presence of two bright galaxies at the center of 
the cluster. The BCG of the galaxy cluster is located at 
a redshift of 0.16, while the second galaxy is at a redshift 



of 0.162. Both galaxies are visible in X-ray, at a projected 
distance of ~ 90 kpc. Bcised on the ratio of their X-ray lu- 
minosities, the smaller galaxy appears to be the BCG of 
a second galaxy cluster with a mass of 10 — 40 percent of 
the mass of MaxBCG J217.95869-f 13.53470, and given the 
similar redshifts of the two galaxies, they are likely interact- 
ing. Using a simple merger model, we estimate the impact 
parameter to be ~ 30 — 100 kpc. 

Using XMM-Newton observations of 

MaxBCG J217.95869-I-13. 53470, we characterize the 
physical properties and the dynamical state of the cluster. 
A single-temperature APEC spectral fit to the -R500 region 
yields a luminosity of ~ 7.9 x lO**^ ergs~^ in the 0.5 — 2 keV 
energy range. We also constructed temperature, electron 
number density, pressure and entropy maps of the cluster. 
To the north and west of the BCG, the temperature 
map shows the presence of hot plasma, with temperatures 
~ 1.5 — 2 times higher than the temperature of its surround- 
ings. The northern hot region lies right to the northwest of 
the smaller galaxy. 

Based on the temperature map, we created projected 
temperature and deprojected electron number density pro- 
files in a sector away from the hot region, and we fitted 
these profiles with a temperature model and a beta-model, 
respectively. From the fit to the deprojected electron number 
density profile we estimate t he mass of the cluster assuming 
a baryon fraction /(, = 0.133 (jConzalez et al.|[2007^ . and find 
a value of 9.2 x 10^^ Mq. Assuming spherical symmetry and 
hydrostatic equilibrium, the mass estimated from the beta- 
model and the temperature model is 5.9 x 10^^ Mq . Exclud- 
ing the central 0.15i?5oo region, we use the scaling relation 
between Ix and M500 to find M500 ~ 1. 1 x lO'^"^ Mp . From 
the _L[0.5 - 2] - Afgoo scaling relation of lPratt et "alT(|2009l) . 
the mass within -R500 is 1.8 x 10^* M©. Considering that we 
are observing a disturbed system, we believe that the mass 
is higher than calculated from hydrostatic equilibrium, so it 
is probably closer to 1 x lO'^'' Mq. 

The regions of high temperature seen in the tempera- 
ture map are also characterized by high entropies. Further- 
more, to the west of the western region, there is a region 
of anomalously high entropy compared to its surroundings. 
Another region of high entropy can be seen to the south- 
west of the BCG. The hot regions, as well as the disturbed 
entropy distribution, support the hypothesis of a merger. 

The southwestern high entropy region coincides with 
the position of the southwestern radio emission visible in 
the GMRT images. The southwestern radio source is con- 
nected to the central radio source by a bridge of faint ra- 
dio emission. Therefore, it is likely that the two sources 
are connected, and the southwestern source traces its ori- 
gin to the center o f the cluster. Based on the radio data of 
Ivan Weeren et al] (|2009l ) and on the evidence of a merger 
from the X-ray observations, we speculate that the south- 
western radio source originated as a bubble from the cen- 
tral AGN, lost significant energy via synchroton and inverse 
Compton emission, and was subsequently revived by the 
merger of MaxBCG J217.95869-I-13. 53470 with the second 
galaxy cluster. This smaller cluster appears to have entered 
the ICM of MaxBCG J217.95869-hl3. 53470 from southeast. 
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the two clusters being currently observed just after the clos- 
est encounter of their BCGs. 

Higher resolution X-ray observations and integral field 
spectroscopy data of MaxBCG J217.95869-f 13.53470 are 
needed to provide more information about the geometry and 
state of the merger. 
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